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Nuclear Forensics and Radiochemistry: Cross Sections 

Robert S. Rundberg 


Abstract: 

The neutron activation of components in a nuclear device can provide useful 
signatures of weapon design or sophistication. This lecture will cover some of the 
basics of neutron reaction cross sections. Nuclear reactor cross sections will also be 
presented to illustrate the complexity of convolving neutron energy spectra with 
nuclear excitation functions to calculate useful effective reactor cross sections. 
Deficiencies in the nuclear database will be discussed along with tools available at 
Los Alamos to provide new neutron cross section data. 
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Radiochemistry: Cross Sections 
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Simple Concept of Cross Section 



V X dt 
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Production of an Isotope by Neutron 

Capture 


The rate of neutron capture is, 

= a(E)<l>(E,t)NA. 

where (j{E) is the capture cross section as a function of neutron energy, (^{E) 
is the neutron flux as a function of neutron energy, and Na is the number of 
target atoms. The number of atoms produced in an experiment is, 

NA+i{E,t)= [ a{E)<P{E, t)NA dt. 

Jo 


For very fast pulses the integrated flux over time, fluence,$(F^) is used and 

Na+i{E) = a{E)^E)NA. 



The Compound Nucleus 


Entrance Compound 
Channel Nucleus 


Exit 

Channels 


A + a —^ (A + a)' 



A + a (Elastic scattering) 


^2 A* + a (Inelastic scattering) 


B+b (Nuclear transformation) 


C+b+c (Nuclear transformation) 


D + d (Nuclear transformation) 


Etc. (Nuclear transformation) 



Neutron Resonances are Levels in the 

Compound Nucleus 



Compound Nucleus 



























Typical Neutron Resonance Spectrum 
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ENDF-B6 

DANCE 



Neutron Energy (eV) 

FIGURE 4. capture cross section measured at DANCE. 


























1/v Cross Section in a Thermal 
Spectrum is Special 


If the neutron velocity distribution is a Maxwell-Boltzmann distribution the flux 
is, 

dSiv) ( m \3/2 / —mv^\ o 

If the neutron reaction cross section has a 1/t? dependence and the flux distri¬ 
bution does not change with time the number of atoms produced in an event 
is, 


Na+i 




vdvdt = aoTioVot. 


Thus the reaction rate is temperature independent when the cross section has 
a 1/v dependence. 







The Breit Wigner Formula 


Neutron cross sections have an inherent l/u dependence when the de Broglie 
wavelength is long compared to the nuclear dimension. This is illustrated in 
the Breit-Wigner formula for neutron capture to a single state in the compound 
nucleus, 


(7(n,j) = TrA^g 


TnT- 


(Er, - Eof + (r/2)" ’ 


where the statistical spin factor g = 


2Zg + l 


(2Z„ + l)(2/t + l) 


, r„ is the neutron width, 


Ty is the gamma width, and T the total width for the compound nucleus, as 


En 0, the Lorentzian becomes a constant IjE^. The deBroglie wavelength, 
A 1/v, and r„ u,thus. 


a(n,'y) l/u. 


The area under a single resonance is. 


A = 


which is proportional to 1/v when ‘C T^y. 







Area Under Lorentzian 


Area under a Lorentzian: details 

To find the area under a Lorentzian (i.e. to integrate from —oo to +oo), make 
the following substitution: 


{u- i/o) = y tan 0 


(Obviously if the Lorentzian is given with different constants, choose 
appropriately: you want to end up with just (1 + idX^O) times some constant in 
the denominator.) Then, 

dv = Y sec^0 = y (1 + tan^ff) d9 


(the latter following from the Pythagoras identity). With the above substitution, 
the new limits on the integral become —nl2 and +^/2 because that's where 
tan 0 goes to —oo and +oo respectively. So now we have 


/ 


y • y(l + ian^ff)d0 
Y^{\ + tan^0) 


which simplifies very nicely to 




= n 



Neutrons per MeV 


Reactor Spectrum 




Flux Distributions 


Watt Spectrum Phys.Rev. 87, 1037(1952) 

P(E) = 0.4865sinh(v^)exp (-£) MeV“' (1) 

See Nuc. Phys. A 957, 289(2017) for a more accurate model of fission neutrons. 
Maxwell Boltzmann Distrbution 

f(E) = 2^(^^y exp i-E/kT) (2) 




Effective Reactor Cross Sections 


In the Wescott convention, the neutron reaction rate for a given nuclide is 
i ?5 = ni'ocfo^ where the effective cross section is given as cfo = <To( 5 f + rs). g is 
a neutron temperature dependent factor that accounts for the deviation of the 
excitation function from a l/n dependence, r is the ratio of epithermal neutrons 
to thermal, s is the ratio of the reduced epithermal integral to thermal, cq is 
the thermal cross sections defined at 2200 m/s neutron velocity. 

s is calculated as _ 

T y/ivao 

where S' is the cross section above the cadium cutoff integrated over sl 1/E 
spectrum with the 1/v cross section subtracted and 

= ( 2 ) 

where / is the ratio of epithermal flux to the thermal flux, fikT is the energy 
below which the epithermal neutron spectrum ends. 







Westcott-g 


The neutron temperature dependent parameter g is defined as the Maxwell 
Boltzmann average reaction rate over the rate for r 1/u absorber, 



( 3 ) 




Activation Products 


Neutron activation products can provide forensic 
information about: 

- Materials used as device components. 

- The location of materials. 

- The neutron spectrum that materials were exposed to. 

- The amount of environmental material swept up into 
the fireball. 

Neutron activation of ubiquitous materials may aid 
in the determination of nuclear yield. 

Neutron capture reactions are sensitive to the 
details of the low energy spectrum. 



Iridium Activation 


Gives information about: 






Neutron energy > 8MeV 

- ''^''Ir(n,2n)''^°lr 
Neutron energy > 1 MeV 

- ''^^Ir(n, n’ 


Low energy neutrons 

- ''^^lr(n,Y)^®2|r 
Fast flux 

- ''^''Ir(n,2n)2 
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Iridium Excitation Functions 


o: 
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Nuclear Decay Characteristics for Ir 

Isotopes 


Isotope 

Produced by 

Half-Life 

Mode of Decay 

Gamma 

Radiations 

(MeV) 

X-Rays 







i89ir 

(n^n)^ 

13.10 d 

EC 

0.245 

Os K and L x-rays 

i90ir 

(n^n) 

11.78 d 

EC 

0.186,0.518,0.605 

Os K and L x-rays 

i92ir 

(n,Y) and (n,2n) 

74.2 d 

P-,EC 

0.308,316,468 

Pt K and L x-rays 
Os K and L x-rays 

193mjj. 

( 1141 ’) 

10.60 d 

IT 


Ir K and L x-rays 

i94ir 

(n,Y) 

19.15 h 

P- 

0.294,0.328 



X-ray Energies (keV) 


Element 

K„ 

Kp 

La 

L, 

Lv 

Os 

61.48> 

71.40> 

8.84 

10.35 

12.09 


62.99 

73.30 

8.91 

10.60 


Ir 

63.28> 

73.55 

9.10 

10.71 

12.51 


64.89 

75.61 

9.17 

10.92 


Pt 

65.11 

75.74 

9.36 

11.07 

12.94 


66.82> 

77.87> 

9.44 

11.25 



Nuclear decay characteristics for and 


Slide 17 























Cross Section (barns) 


Gallium Cross Section Needs 
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10-10 

-T 


10-5 

-T 


1 

T 


10 ^ 













Cross Soction (barns) 


Another Gallium Cross Section Need 
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Cross Section (barns) 


Mistakes Have Been Made 

























Ga-68 Decay Scheme 


**Ga c Decay l999BeZQ,l999BeZS 


Decay Scheoie 


Intensities: Ky^ce) per 100 parent 
decays 



67.71 min 

/ S1<5*17 

%e+%B+= 100 

<?‘=2921.1*’ 


Ee 

IP* 

le 

bog f 



0.0104 

5.113 



0.0955 

5.754 



0.230 

5.87 9 


0.000258 

0.0328 

6.898 


1.19 

1.80 

5.489 

stable 2912 

87.72 

8.92 

5.1942 


































Annihilation Gamma Ray Summing 



1 + Ot-] 


( 1-2 


4^2 

1 + 0t2 



Summing with annihilation radiation is a 
special case because the photons are spa¬ 
tially correlated. The probability of detect¬ 
ing a 511 keV photon is doubled. 




Device for Advanced Neutron Capture 

Experiments, DANCE 




DANCE 




Layout of the DANCE Time-of-Flight 

Spectrometer 
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DANCE Location in the Lujan Center 
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Properties of BaF2 Scintillators 



PMT Ajjode — 


The total light output from BaF 2 arises 
from two scintillation mechanisms, 
which are differentiated by their 
respective decay times (e.//. fast and 
slow). 

BaF 2 suffers from an intrinsic a 
activity (averaging ~225 cnt/s/l for 
crystals in the DANCE array) resulting 
from unavoidable Ra contamination (a 
chemical homologue of Ba). 



Digitizer Channel (2 ns / channel) 















Scintillation Properties 


Scintillation Properties 

Fast 

Slow 

Emission Maximum A (nm) 

220 

310 

Primary Decay Time (ns) 

O.G - 0.8 

G30 

Light Yield (photons/MeV-g) 

2.5 X 10'* 

G.5 X 10'* 

Photoelec. Yield (% of Nal(TI)) 

5 

IG 






DANCE Pulse Shape Discrimination 


Relative intensities of fast and slow 
components exhibits a strong 
dependence on the type of radiation 
detected. 



a-particles deposit energy almost 
exclusively through the slow component 
mechanism. Comparison of the fast and 
slow component integrals provides a 
means to reject the background (\ 
activity. 



r^sum (i‘'h- units) 


Digitizer Channel 




















Neutrons/(cm^ eV Tq) 


DANCE Neutron Flux 



10 ^ 10 ® 10 '* 10 ® 
Neutron Energy (eV) 







Counts 


DANCE Backgrounds 



Summed gamma-ray energy spectra for reactions with neutron energies between 1 and 10 keV and gamma-ray multiplicity > 2. 








U-238 n, gamma Resonances 



Energy (eV) 



















Cross Section (barns) 


U-238 n,gannnna 
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Moxon 06 

A 

Adamchuk 77 

O 

Fricke 71 
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Linenberger 44 
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10 ^ 


Energy (eV) 


10 


Measured ‘^*U(n, 7 ) cross section compared to previous measurements by Moxon [3 1 Adamchuk 32], Fricke [30] and Linenberger 33]. 



























Cross Section (barns) 


U-238 Unresolved Region 



Energy (eV) 


Measured ^^U(n, 7 ) cross section compared to previous measurements by Kazakov [3vj, Poenitz [35], Moxon [3i i, Adamchuk [32], Fricke [30], and 

















Counts 


U-238 Gamma Ray Spectrum on a 

Single Resonance 



Spectrum of 7 rays with Mcj = 2 from the decay of the 36 eV resonance in ^^^U(n, 7 ) compared to the blank background subject to the same data 
reduction gates. 


















Comparison to DICEBOX 



Spectrum of 7 rays from the decay of the four resonances in ^^U(7i,7), Md = 2 , compared with dicebox calcuiations (GLO) with parameters as 
described in the text. 







□ICEBOX Multiplicity 3 
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Spectrum of 7 rays from the decay of resonances in ^^^U(n, 7 ), = 3, compared dicebox calculations. 









U-238 Study Conclusion 


The gamma ray spectra observed required Ml 
transitions in the DICEBOX cascades. 

The Ml transitions are due to scissor-mode 
vibrations as opposed to E2 rotations 



ENDFB-VI Problem Solved 
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FIGURE 5. capture cross section measured at DANCE, The data is normalized to the thermal 

neutron capture cross section. 

























